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Abstract 

We reanalyze the magnetization data collected on E^S^Cai-zYzC^Og.^ sam- 
ples (Kim at al, Phys. Rev. B 72, 64525 (2005)) and argue that the method, which 
was used for the analysis of equilibrium magnetization data, is not adequate to 
the experimental situation. As a result, the temperature dependencies of the upper 
critical field H&iT) and the magnetic field penetration depth A(T), obtained in 
this work, are misinterpreted. Using a different approach to analysis, we demon- 
strate that the normalized H&iT) curves are rather different from those presented 
in the original publication and do not follow predictions of the Werthamer-Helfand- 
Hohenberg theory. Another important observation is that the H&iT) dependencies 
for two samples with different levels of doping are qualitatively different. 
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1 Intrtoduction 



In this paper, we reanalyze mixed-state magnetization M(H, T) data that were 
collected on several polycrystalline Bi 2 Sr 2 Ca 1 _, r Y, r Cu20 8+( 5 (Bi-2212) samples 
and presented in Ref. pQ. The main reason is that these data contain some 
hidden information, which can be discovered by an appropriate analysis. We 
also argue that the method, which was used for the analysis of equilibrium 
magnetization data in p], is not adequate to the experimental situation and 
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the doping dependence of the zero-temperature value of the upper critical field 
H c2 (0), presented in pQ, should be considered as unjustified. 

As was demonstrated in [215PI] . all numerous high-T c superconductors (HTSC) 
may be divided into two groups. The dependencies of the normalized upper 
critical field h c2 on T/T c for HTSC's belonging to the same group are prac- 
tically identical, while they are distinctly different between the groups. Ac- 
knowledging the fact that the larger group includes a huge variety of different 
HTSC compounds, while the second one is rather small, we shall denote the 
corresponding h c2 (T/T c ) curves as typical and unusual, respectively. Quite sur- 
prisingly, the results for two Bi-2212 samples with different levels of doping, 
investigated in [1] and analyzed in this work, perfectly match the correspond- 
ing h c2 (T/T c ) curves for the two above mentioned groups of HTSC's. As well 
as we are aware, it is the first observation of such a behavior in Bi-based 
HTSC's. 

There are several theoretical approaches, which are usually employed for eval- 
uation of H c2 from experimental magnetization data [5.6.7J. All these models 
assume conventional superconductivity (an isotropic superconducting order 
parameter) and a uniform sample with a zero demagnetizing factor. Neither 
of these conditions is satisfied in polycrystalline HTSC's. Because the differ- 
ences between theoretical assumptions and experimental situations are rarely 
discussed in the literature, we consider them in some details. 

(i) Demagnetizing factor. If the sample magnetization M is much smaller than 
an applied magnetic field H, demagnetizing effects are usually neglected. This 
is not correct. If a demagnetizing factor n ^ 0, the sample magnetization can 
be written as M nj t Q = (1 — n)M n=0 (in the case of 4ttM <C H), i.e., M — > 
for n — > 1, independent of an applied magnetic field, temperature or the 
nature of the sample. Demagnetizing effects are not important if only relative 
variations of the sample magnetization are considered. In Ref . [TJ , however, the 
Hao-Clem [5] and the vortex fluctuation [6fT] models were employed for the 
analysis of experimental data. The absolute values of M enter both of these 
models and neglecting demagnetizing effects may result in misinterpretation 
of experimental results. 

(ii) Pairing symmetry. Symmetry of the order parameter is also important. In 
the case of unconventional <i-pairing, which is expected in HTSC's, the dis- 
tribution of the order parameter around vortex cores and the corresponding 
contribution to the free energy is different from that for conventional supercon- 
ductors. This is why theoretical calculations based on conventional s-pairing 
should be used with caution if they are applied for the analysis of experimental 
data collected on unconventional superconductors 

(iii) Polycrystalline samples. HTSC's are highly anisotropic. In such materials, 
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if the direction of an external magnetic field does not coincide with one of the 
principle axes of the crystal, the magnetic induction in the sample is not ex- 
actly parallel to the applied magnetic field. In samples consisting of randomly 
oriented grains, this leads to an additional free energy and may influence the 
sample magnetization. It should also be noted that, because magnetizations 
of different grains are different, there is some magnetic interaction between 
the neighboring grains. The situation is even more complex at higher tem- 
peratures. Indeed, according to calculations of Brandt [8J, the magnetic field 
dependence of M is a linear function of In if (London limit) only in magnetic 
fields H < 0.1H c2 (see also Fig. 3 in Ref. p|). At temperatures, T > 0.8T C , 
the upper limit of the magnetic field range is usually higher than this value. 
In this case, deviations of M(H) from the predictions of the London model 
have to be accounted for and a simple averaging, which was proposed in Ref. 
[7j and used in Ref. [I], is not applicable. 

The model of thermal fluctuations of vortices [6], which was used for the anal- 
ysis of experimental data in Ref pQ, is based on an experimental observation 
that, in the case of layered HTSC compounds, there is a temperature T* < T c , 
at which the sample magnetization does not depend on an applied magnetic 
field [10J. As well as we are aware, the model [Bj is the only theory describ- 
ing this feature. Using this approach, one may evaluate the magnetic field 
penetration depth A at T = T*. The only parameter s, entering the expres- 
sion for A(T*), represents the distance between superconducting layers. This 
parameter may be independently evaluated as s = — (^T*)/(f M*) is 
the Boltzmann constant, and $o is the magnetic flux quantum). However, the 
ratio —T*/M* is always smaller than the theoretically predicted value and, 
contrary to the theory, T*/M* is practically independent of s [11] . This is why 
the model j6] may be considered only as a qualitative approach to the problem 
and the resulting value of A(T*) may be different from the actual magnetic 
field penetration depth. 

One of possible reasons of the above mentioned disagreement between the the- 
ory and experiments is that thermal fluctuations of vortices, considered in [6], 
is not the only fluctuation effect that may contribute to the sample magneti- 
zation. As was discussed in [2] , fluctuations of T c throughout the sample vol- 
ume, which cannot be avoided in HTSC's, should also play an important role. 
Cuprates are non-stoichiometric, which makes them intrinsically inhomoge- 
neous materials. Fluctuations of chemical composition cannot be smaller than 
the corresponding statistical fluctuations in distribution of non-stoichiometric 
components. In real samples, however, chemical fluctuations are even stronger 
than the corresponding statistical numbers. Taking into account that T c is 
strongly dependent on the level of oxygen or other dopants, we may conclude 
that the critical temperature must be spatially dependent. 

At T above the bulk critical temperature, our sample can be considered as 
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superconducting inclusions (grains) imbedded in normal metal. The same sit- 
uation can be observed at T < T c in magnetic fields H > H c2 . At temperatures 
well below T c and in H -C H c2 , the above mentioned non- uniformity of the 
sample is not important and it may only lead to weak variations of the vortex 
density. In this case, the sample magnetization will correspond to an averaged 
value of H c2 . However, closer to the H c2 (T) line, the situation is different. 
Indeed, if H > H c2 (T), we have only superconducting grains and the diamag- 
netic moment of the sample should be proportional to H (see Eq. (4.12 in 
[T2]). In magnetic fields close but below H c2 (T), the superconducting order 
parameter A is small in the bulk of the sample (A — > for H —>■ H c2 ), while it 
is substantially higher in the regions where local T c is higher. Considering the 
sample magnetization, we may conclude that it consists of two diamagnetic 
contributions: M gr (from inclusions with higher T c ) and M ms (from the mixed 
state between the inclusions). Both these contributions are approximately lin- 
ear on At the same time, the derivatives dM gr /dH and dM ms /dH have 
opposite signs. At low temperatures, M ms 3> M gr . However, M ms vanishes 
at the bulk value of T c , while M gr remains non-zero up to somewhat higher 
temperatures. In other words, it must be a temperature T* < T c , at which 
(dM ms /dH+dM gr /dH) = 0, i.e., M = M ms +M gr is temperature independent. 

As we could see, spatial variations of the superconducting critical tempera- 
ture in inhomogeneous samples result in an effect similar to that of thermal 
fluctuations of vortices considered in [B]. It is quite likely that both these ef- 
fects are important and this is why the interlayer distance s, calculated from 
experimental values of T* and M* as s = — (A;bT*)/($oM*), is in disagree- 
ment with its experimental results. It is also possible that in polycrystalline 
samples, which are expected to be less homogeneous than single crystals, T* 
and M* are entirely determined by spatial fluctuations. 



2 Model 



Here, we use a completely different method of analyzing of magnetization 
data. In this scaling approach, developed in Ref. [2], no particular M(H) 
dependence is assumed a priori, and it can be applied to single crystals as 
well as to polycrystalline samples, independent of the pairing mechanism or 
the sample geometry j3]. The disadvantage of this analysis is that it does not 
provide the absolute values of H c2 but only its relative temperature variations. 
This is the prize to pay for its universality. However, if the value of H c2 at any 
temperature is established, a whole H c2 {T) curve is obtained. 

The scaling procedure is based on the assumption that the Ginzburg-Landau 

1 According to [8], M ms (H) is approximately linear function down to H 0.3H c2 . 
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parameter k is temperature independent. In this case, the mixed-state mag- 
netic susceptibility may be written as 

X(H,T) = X (H/H c2 ), (1) 

i.e., the temperature dependence of % is om Y due to temperature variation of 
H c2 . Eq. (1) is already sufficient to establish a relation between magnetizations 
at two different temperatures [2J 

M(H/h c2 ,T )=M(H,T)/h c2 , (2) 

where h c2 {T) = H c2 (T) / H c2 (T ) is the upper critical field normalized by its 
value at some arbitrary chosen temperature T < T c . This equation is valid if 
the diamagnetic response of the mixed state is the only significant contribution 
to the sample magnetization. Considering HTSC's, however, we have to take 
into account their noticeable paramagnetic susceptibility Xn in the normal 
state and its dependence on temperature. In order to account for Xn(T), we 
have to introduce an additional c (T)H term in Eq. (2). According to [2J, the 
resulting equation connecting M(H,T ) and M(H,T) may be written as 

M(H/h c2 , T ) = M(H, T)/h c2 + c (T)H (3) 

with 

c (T) = X n(T) - Xn(T ) (4) 

We note that Eqs. (1) and (2) are rather general and they can be obtained 
from any model based on the Ginzburg-Landau theory, including the so-called 
nonlocal London theory and the Hao-Clem model. At the same time, as was 
discussed in [13111"] . these relations remain valid even if M(H) is different from 
predictions of the conventional Ginzburg-Landau theory and, therefore, they 
are applicable to unconventional superconductors as well. 

It should be noted that at temperatures close to T c , some additional contribu- 
tion to M arises from fluctuation effects. As was discussed in [2J, the second 
term in Eq. (3) may also account for this contribution. However, because this 
contribution is not exactly proportional to H, it can be accounted for only 
partially. 

Eq. (3) can be used as the basis for the scaling procedure. The adjustable pa- 
rameters h c2 (T) and co(T) are obtained from the condition that M(H/h c2 , T ), 
calculated from data collected at different temperatures, collapse onto a single 
master curve, which represents the equilibrium magnetization at T = T [2J. 
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As a result, the temperature dependence of the normalized upper critical field 
h c2 (T) is obtained. Extrapolating the resulting h c2 (T) curve to h c2 = 0, we 
can also obtain the value of the zero-field critical temperature. 

It was demonstrated that this scaling procedure works quite well and may be 
used in order to reliably obtain temperature dependencies of the normalized 
upper critical field h c2 from equilibrium magnetizations measured at different 
temperatures [2|3|9|13f 14lll5f 16] . In the following, we use M e ff(H) to denote 
M(H/h c2 ,T ) calculated using Eq. (3) in order to distinguish it from directly 
measured magnetization data. 



3 Analysis of magnetization data 

Figs. 1(a) and 1(b) show scaled magnetization curves for two Bi-2212 samples 
with different Yttrium contents. As may be seen, the quality of scaling is 
almost perfect in both cases and deviations of individual data points do not 
exceed the accuracy of the original data as they can be taken from the figures 
presented in pQ. 

The resulting temperature dependencies of the normalized upper critical field 
for these samples are shown in Figs. 2 and 3. In order to demonstrate a weak 
interference between the two fit-parameters, we repeated the scaling proce- 
dure assuming cq = 0. As may be seen in Fig. 2, the difference between the 
two sets of data-points is insignificant (see also [15]). At the same time, be- 
cause the normal-state paramagnetism in HTSC's exists and it is temperature 
dependent, we do not see any reason to neglect c (T). 

The results for a sample with the Y-concentration x = 0.2 are shown in Fig. 
2. The solid line represents the "typical" h c2 (T/T c ) dependence, obtained in 
Ref. [2J, and fitted to the data points by adjusting T c . The value of T c = 92.6 
K, evaluated in such a way, is close to the value given in the original pub- 
lication!!] As may be seen in Fig. 2, the data points follow the solid line 
quite closely clearly demonstrating that this sample belongs to the above men- 
tioned larger group of HTSC's. We also note that because the h c2 {T) curve 
for this group of HTSC's is quite different from predictions of the Werthamer- 
Helfand-Hohenberg (WHH) theory [T7| (the dashed line in Fig. 2), the zero- 
temperature upper critical field H c2 (0), evaluated by employing this theory 
will be well above its real value. 

The h c2 (T/T c ) curve for a sample with x = 0.1, as shown in Fig. 3, is quite dif- 



2 Because the zero-field superconducting transition in HTSC's is rather broad, the 
value of T c , evaluated as it was done in [1], may be considered only as approximate. 
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Fig. 1. The scaling results for two Bi-2212 samples with different doping levels. 
Magnetization data are taken from Figs. 1 and 2 of Ref. [TJ. 

ferent and practically coincides with those for several other HTSC compounds 
belonging to the smaller group of HTSC's. Our evaluation of T c = 84.2 K for 
this sample is practically the same as the value provided in Ref. [T]. We note 
that, although in this case h C 2(T/T c ) is closer to the WHH theory, the differ- 
ences are still significant and this theory should not be used for the evaluation 
of H c2 {0). 



4 Discussion 



The observation that, considering temperature dependencies of H c2} all nu- 
merous HTSC's may be divided into two groups is remarkable. However, it 
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Fig. 2. The normalized upper critical field h c2 as a function of temperature for a 
sample with x = 0.2. For comparison, we also show h C 2(T) obtained if the scaling is 
based on Eq. (2) (co(T) = 0). The solid line represents the best fit of the "typical" 
h C 2(T/T c ) curve, obtained in Ref. [;2J, to the data points (see text for details). The 
dashed line represents to the WHH theory. 




Fig. 3. The normalized upper critical field h C 2 as a function of T/T c for a sample 
with x = 0.1. The h C 2(T/T c ) curves for several samples, belonging to the smaller 
group of HTSC's |3|13j . are shown for comparison. The solid and the dashed lines 
are the same as in Fig. 2. 

is even more surprising that no any intermediate h c2 {T /T c ) was observed so 
far. Because both typical and unusual h c2 (T/T c ) curves were observed in the 
same families of HTSC's [2f3]P| , one may assume that the level of doping is 
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essential. A similar conclusion may also be drown from the results presented 
in this work. This is why it would be extremely interesting to study the tran- 
sition from one type of the h c2 (T/T c ) dependence to the other systematically. 
We also note that in this Y-doped Bi-2212, an overdoped sample fall into 
the smaller group (unusual h c %{T/T^} curves), while in the case of Y-123 only 
underdoped samples behave in such a way (see Fig. 3). 

Finally, we briefly discuss the results of Ref. [1]. The main result is the de- 
pendence of the zero-temperature upper critical field H c2 (0) on the doping 
level. In order to obtained this plot, the authors had to extrapolated their 
H c2 (T) data to T = 0. However, as we argue below, neither the H c2 {T) curves 
nor their extrapolations can be considered as reliable. We also note that the 
H c2 (T) curves presented in [Tj are quite different from our results (see Figs. 2 
and 3). 

As may be seen in Fig. 4 of Ref. pQ, the H C 2(T) data points cover a rather 
narrow range of T/T c values. In such cases, employing theoretical expressions 
for extrapolation of experimental data is justified only if it is esyablished that 
the corresponding theory is quantitatively applicable. As well as we are aware, 
there are no experimental proofs that the modification of the WHH theory, 
proposed in [18J, can be used for description of HTSC's. Moreover, there are 
no theoretical reasons to expect this. Indeed, the theories [TT1[T8] are based on 
the conventional BCS theory and their applicability to unconventional super- 
conductors is questionable. We also note that, considering the H c2 (T) data 
presented in Fig. 4 of Ref. [T], one can easily see that the data-points are in 
disagreement with the theoretical curves, which were used for their extrapo- 
lation. 

Another warning is that the presented H c2 (T) data do not show a tendency 
to vanish at T = T c . This is a strong indication on some drawbacks in the 
theory, which was used for evaluation of H c2 from magnetization data. As we 
discussed in the Introduction, this theory neglects spatial fluctuations of T c 
and this is why it is quite likely that it is not applicable to real samples. In 
this case, the data, presented in Fig.4 of Ref. [1], do not represent H c2 and 
the conclusions made in this work are not actually based on the presented 
experimental data. 

In conclusion, using an alternative approach to the analysis of experimental 
data presented in Ref. [I], we demonstrate that, depending on the Y-content, 
Bi 2 Sr 2 Cai_ :r Y 2; Cu 2 8+( 5 samples may have qualitatively different temperature 
dependencies of the upper critical field. This result is in agreement with previ- 
ous observations of similar behavior in other HTSC compounds [31TT3] . Another 
important point is that our H c2 (T) curves are rather different from the results 
of the original publication. We argue that this disagreement between two anal- 
ysis of the same experimental data is due to non-applicability of the vortex 
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fluctuation model to description of real HTSC samples. 

This work is partly supported by the Swiss National Science Foundation. 
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